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We study the proximity coherence in a mesoscopic normal-metal film (N) in contact with a 
superconductor (S). Accounting for a repulsive interaction between the electrons in the normal 
metal, we find an enhanced local density of states close to the NS interface. The sharp peak in 
the density is pinned to the Fermi energy and leads to spontaneous paramagnetic interface currents. 
The induced orbital magnetic moments exhibit the characteristic features of paramagnetic reentrance 
observed in normal-metal coated superconducting cylinders [Phys. Rev. Lett. 65, 1514 (1990)]. 



A normal metal in contact with a superconductor ex- 
hibits the phenomenon of proximity — the superconduc- 
tor exports its coherent state across the interface into 
the normal metal. On a microscopic level, this phe- 
nomenon is described through the Andreev reflection of 
the normal-metal quasi-particles at the normal-metal 
superconductor (NS) interface, converting normal- to su- 
percurrent. Proximity superconductivity exhibits a rich 
phenomenology and has attracted considerable interest 
recently [Q. A particularly puzzling finding is the ultra- 
low-temperature reentrance observed in normal-metal 
coated superconducting cylinders where, contrary to 
expectation, the fully diamagnetic cylinder develops a 
paramagnetic response at low temperatures. Recently, it 
has been speculated that a novel kind of persistent cur- 
rent states circling the cylinder might be responsible for 
this phenomenon Q , but closer inspection of the experi- 
mentally measurable quantities reveals that the predicted 
effect is too small Q. In this Letter, we demonstrate that 
the presence of a repulsive electron-electron interaction 
in the normal metal naturally leads to the appearance of 
a paramagnetic instability at very low temperature, of- 
fering a possible explanation of the reentrance effect in 
the NS cylinders. 

To be specific, we shall consider a clean normal-metal 
slab of thickness d (0 < x < d), in perfect contact with a 
bulk, conventional superconductor. The proximity ef- 
fect is mediated by the Andreev reflection at the in- 
terface with the superconductor, which transforms inci- 
dent electrons into back-reflected holes, thus binding the 
quasi-particles states to the normal layer for E < A5. 
In the usual free electron gas description of the normal 
metal, the Andreev bound states are found at E n — 
hv x (2n + l)7r/4d (n = 0,1,..; v x = VFCOsd) producing 
a linear suppression of the DOS M N(E) ~ N Ed/hv F 
close to the Fermi level E = (N = mkp /h 2 ir 2 ). In 
the following we assume that the electron-electron in- 
teraction in the normal layer, which follows from the 
delicate balance between the phonon-mediated- and the 
Coulomb-interaction, is repulsive. As a consequence, a 



finite order parameter A(x) is induced in the metal, op- 
posite in sign as compared to As in the superconduc- 
tor, see Ref. 0. The NS junctions then behaves like a 
Josephson junction with a phase difference 7r, trapping 
quasi-particle states at the Fermi energy close to the NS 
interface. The local density of states N(E,x) exhibits 
a peak at zero energy on top of the Andreev density of 
states, as shown in Fig. [l]. This peak involves a macro- 
scopic number of states with density n p ~ k 2 F / d, which 
in the following we call the 7r-states. 

The change in the DOS crucially affects the response 
of the proximity metal. The linear current response j [A] 
can be divided into two contributions j — jdia~t~Jpara; 

the 

diamagnetic current jdia = — (e 2 n/mc) A giving the rigid 
response of the bulk density n = k F /2>n 2 and the para- 
magnetic current j para following from the deformation of 
the wavefunction at the Fermi surface , 

e 2 n . f / df\ N(E) 

for slowly varying fields A (/ is the Fermi occupation 
number). While in a bulk superconductor the param- 
agnetic current is quenched by the energy gap at low 
temperatures producing a net diamagnetic response, the 
paramagnetic current of a bulk normal metal cancels 
the diamagnetic current exactly. In the non-interacting 
metal under proximity, the linear density of states sup- 
pression N(E) oc E is still sufficient to suppress the para- 
magnetic current at zero temperature J^]. Including a 
repulsive interaction places the system in the opposite 
limit: The sharp DOS peak at the Fermi level produces a 
paramagnetic signal which cwercompensates the diamag- 
netic response. Such a paramagnetic response naturally 
leads to an instability: The free energy SF = —cjSA < 
can be lowered via a non-zero magnetic induction induced 
by spontaneous currents along the NS interface. The in- 
terface currents are associated with an orbital magneti- 
zation M(T) producing a low-temperature reentrance in 
the magnetic susceptibility. 

In the following we present a quantitative analysis 
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of the paramagnetic instability induced by the 7r-states. 
The magnetic induction B z {x) parallel to the surface is 
described by the vector potential A y (x) which drives the 
currents j v (x). The electron-electron interaction in the 
superconductor is accounted for by an effective coupling 
constant Vs < and similarly Vn > in the normal 
metal, see also Refs. f) 10 1. Two self-consistency prob- 
lems have to be solved: First, we evaluate the order pa- 
rameter A(x) accounting for the different coupling con- 
stants in the superconductor and the normal metal, and 
obtain the local DOS N(E,x). Second, we determine 
the current functional j [A] which we solve together with 
Maxwell's equation to find the spontaneous interface cur- 
rents. 
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FIG. 1. Local DOS N(E,x) at the NS interface x = and 
at the metal boundary x — d (Ea ~ Tivf/cL), as it follows 
from the self-consistent solution of the Eilenberger equation, 
Eqs. (Q) and (|^), for a thickness d = IOHvf / As and the 
coupling constants Vs = —0.3 and Vn = 0.1. Inset: Spa- 
tial dependence of the order parameter A(x) and local DOS 
N(E = 0, x) at the peak energy (S: x < 0, N: x > 0). 

We use the quasi-classical description following from 
the Eilenberger equation for g, 

- hv x d x g = [{hu + iev y A y (x)/c}T 3 + A(x)ri,g] , (2) 

where the 2x2 matrix g contains the Green's functions 
gu{x, v x ) and fui{x,v x ) (to Matsubara frequency, v x = 
vp cosd, fj Pauli matrices, see Ref. |ll]]). Eq. (||) is com- 
pleted by the self-consistency relation for the pair poten- 
tial ((..) is the angular average), 



A(x) 



-VN 7tTJ2(L(x,v x )). 

w>0 



(3) 



The self-consistent numerical solution of Eqs. (g) and 
(||) is shown in the inset of Fig. |l]. The course of the 
order parameter in the normal layer is asymptotically 
given by A(x) ~ —VnNqTivf/x, as expected from the /- 
function in the non- interacting case Vn = 0. A(x) decays 



from a value ~ — \Vn/Vs\As at the NS interface, to ~ 
— VmNqTivf/cL at the outer boundary. Close to the NS in- 
terface, the local DOS N(E,x) = N Re{(g^iE+s(x, v x ))} 
exhibits a pronounced peak at zero energy, as shown in 
Fig. [l]. In order to proceed with analytical results, we 
approximate the order parameter by a step function, 



A(x) = 



A, 
-A 



N- 



x < 0, 

< x < d, 



where A^r oc Vjv enters as a parameter. The Green's 
function in the normal layer x > can be determined 
exactly and takes the form 



g u {x,v x 



hui sinh {x{d) — 7] + An cosh [x(d — x)] 
TiQ cosh lx(d) — 7] 



(4) 



where x( x ) — 2Clx/v x , h Q = A N + h to , and tanh7 = 
Ajv /M7 (we consider the limit A5 3> Ajv, hu). The sec- 
ond term in (Q) describes the 7r-states at the NS interface. 
The poles of the Green's function at Tllo = —iE + yield 
the bound state energies. While for E ;§> An the An- 
dreev states of the free electron gas are down-shifted by 
SE n w -2AAr/(2n + 1)tt (n = 0,1,..), below the gap 
E < An we find the 7r-states at 
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exponentially close to Fermi energy. All trajectories with 
v x = vp cosz? <c ANd/h possess a bound state at E w 0, 
thus producing the macroscopic weight of the zero energy 
DOS peak: For A^ > Tivp/d the number of 7r-states per 
unit surface iV sur f is equal to the number of transverse 
levels iV sur f ~ k F , while for A^ < hvp/d it is reduced to 
IVsurf ~ kp(ANd/TiVF) 2 via the reduction of the available 
solid angle cos?? < ANd/hvF- 

We derive the current-field relations at low tempera- 
tures, assuming fc^T <C Tivf /d, An- This implies a ther- 
mal length £,n(T) = hvp/^nkBT larger than the thick- 
ness d and no thermal smearing on the scale An- Only 
the trajectories with cosi? < ANd/hvF contribute to the 
current at low temperatures. We describe them in the 
limit hv x /ANd — > by 
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The current in the presence of a slowly varying vector 
potential A, follows from Eq. (||) after replacing uo by 
w + ieVyA/hc and inserting it into the current expres- 
sion j y (x) = ieN 2irTJ2u;>o( v y9( x , v x,Vy)). In addition 
to the usual diamagnetic current jdia = —(c/4ir\ 2 ) A 
[A = (mc 2 /47rne 2 ) -1 / 2 denotes the London length], we 
obtain the paramagnetic current 
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in the limit evpA/c, k B T <C Ajv. Here, $o = Trfi-c/e de- 
notes the flux quantum, = ^vp/2A^ gives the extent 
of the 7r-states, and under the assumption Ajy > hvp/d 
we have a = 1. At temperatures fcsT ^> evpA/c the 
paramagnetic current j pa ra ~ (c/ \ 2 )(An /}cbT)A is lin- 
ear in A and oc 1 /T, a signature of the thermally smeared 
zero energy DOS peak, and competes with the diamag- 
netic current on the scale fjy. At T — > 0, Eq. (]?]) is non- 
linear in the field and generates the spontaneous para- 
magnetic current. This paramagnetic interface current 
results from the energy splitting of the 7r-states in the 
field, E m ±evpA/c, allowing the system to gain en- 
ergy by shifting the DOS below the Fermi surface. For 
Ajy < hvp/d, the paramagnetic current is reduced by 
the factor a = (A^d/hvp) < 1 in Eq. (^). The surface 
current I = J jdx ~ a 2 c$o/A 2 is in agreement with the 
current estimate I„ ~ N sur {evp based on the number of 
7r-states at zero energy. Eq. (0) thus always produces a 
net paramagnetic response at low temperature and fields. 
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FIG. 2. Magnetization M(T, H) curve at various temper- 
atures (in units of Ajv/fce), f° r A = 0.3d and — d. The 
two meta-stable branches P± exhibit a spontaneous magne- 
tization in zero field, the diamagnetic branch D is unstable. 
Inset: Zero field magnetization Mo(T). 

The evaluation of the induced magnetization re- 
quires the self-consistent solution of Maxwell's equation 
— d 2 A(x) = A-Kj{x)/c together with the current func- 
tional j[A(x)]. The solution of the screening problem re- 
quires the full dispersive relation between j(q) and A(q), 
which in the proximity effect typically features a non- 
local current-field dependence on the scale £jy fTl] , Eq. 
(^) represents the long wavelength limit q — ► 0. For sim- 
plicity, here we use Eq. (^) under the assumption of local 
response (accounting for the full non-locality does not af- 



fect the qualitative nature of the final results, as we have 
checked numerically). 

The magnetization curves M(T, H), which follow from 
Eq. (Fjj) are shown in Fig. |^. Approaching from large 
fields, Fig. H shows two paramagnetic branches P± with 
a linear diamagnetic slope exhibiting a spontaneous mag- 
netization in zero field. They result from the superposi- 
tion of the paramagnetic magnetization Mq(T) and the 
Meissner response to the applied field H. As the field is 
decreased (increased) past H — 0, the branch P + (P_) 
becomes meta-stable. The spontaneous magnetization 
M (T) appears below a second order transition point 
and saturates at low temperatures, as shown in the in- 
set of Fig. 0. The magnetization curve includes a dia- 
magnetic branch D, which arises from the competition 
between the paramagnetic instability and the thermal 
smearing and is thermodynamically unstable. 

In the following we give a semi-quantitative analysis 
of the magnetization M = J dx M(x)/d, first at zero 
temperature and field [Mo], proceeding to finite temper- 
atures [M (T)], and finally including an applied mag- 
netic field H [M(T, H)]. The boundary conditions are 
given by A(0) = and d x A(d) — H. We concen- 
trate on the most relevant limit where £%,d 3> A. At 
T = 0, according to Eq. (f7j), the paramagnetic interface 
current j ~ ac$o/A 2 £^ remains unscreened until being 
matched by jdia ~ — cA/X 2 , producing a vector potential 
A ~ a$o /£jy on the scale A. The vector potential A sat- 
urates beyond A, as the para- and diamagnetic currents 
cancel each other. Assuming that the 7r-states extend up 
to the outer metal surface (a < 1), the induced magne- 
tization M = A(d) /And is given by 
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We note that although the spontaneous currents increase 
as Ajv > hvp/d (a = 1) they are screened exponentially 
beyond the extent of the 7r-states in this limit, giving 
a magnetization M ~ (& /£%d) exp[— (d — £%)/X\. We 
assume a < 1 in the following. 

At finite temperature, the spontaneous magnetization 
is suppressed by the factor a,r cta,n(ev f A/ cnlisT), which 
itself depends on the magnetization via A ~ Md, imply- 
ing the implicit equation 
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The spontaneous magnetization appears below a second 
order transition at fcgT*' ~ a An, saturating at low tem- 
peratures, as shown in the inset of Fig. ^|. The transition 
temperature is equal in magnitude to the energy splitting 
of the DOS peak E ~ evpA/c ~ aA^. 

Under an applied magnetic field H , the Meissner cur- 
rent jdia screens both the spontaneous interface current 
and the applied field. At zero temperature we deal with 
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a linear problem and the magnetization is given by the 
superposition M(H) = M +xH of the spontaneous mag- 
netic moment Mo and the Meissner response \H ■ As the 
temperature increases, Mq(T) decreases and the meta- 
stable regime shrinks. At T > T C M the spontaneous mag- 
netization in zero field has disappeared, the signature 
of the paramagnetic currents remains, however, reduc- 
ing the diamagnetic susceptibility \ & t small fields. At 
large temperature T^T C M we recover the pure Meissner 
response. 

Note that the two meta-stable branches P+ and P_ in 
the magnetization curve, see Fig. [^, imply a first order 
transition with changing field at H = 0. The first order 
transition is similar to the magnetic breakdown occurring 
in the same system at large fields between the fully dia- 
magnetic phase and a field penetration phase [fl2|| . The 
rotation of the magnetic moments to the energetically 
more favorable polarization will show the hysteretic be- 
havior typical for a first order transition. The transi- 
tion from P + to P_ implies a paramagnetic slope in the 
thermodynamic dc-magnetization (M(T,H)), which will 
link the meta-stable solutions P± in Fig. || and cross the 
origin at M(H = Q) = 0. In summary we find that 
on approaching T C M from above, the diamagnetic sus- 
ceptibility Xdc — (M(T,H))/H is reduced, exhibiting a 
low-temperature reentrance. Below T C M , the spontaneous 
interface currents produce a net paramagnetic suscepti- 
bility Xdc- 

In the following we discuss our results in the context of 
the experiments by Mota and co-workers, who have mea- 
sured the magnetic response of normal-metal coated su- 
perconducting cylinders at low temperatures [0,|l3j- Re- 
cent studies have established a quantitative understand- 
ing of the magnetic response of these samples at higher 
temperatures, explaining the sensitivity of the screen- 
ing to small impurity concentration [ [1 l|Jl 4f| due to the 
nonlocality and the magnetic breakdown at finite field 
|i2|| . The Nb-Ag and Nb-Cu cylinders show an anoma- 
lous paramagnetic signal in the magnetic response in the 
low-tempcraturc - low-field corner of the H — T phase 
diagram PJl^|. A direct comparison with our theory 
requires the magnetization curve (M(T, H)) which has 
not yet been measured. The observed dc-susceptibility 
Xdc(T) as a function of temperature shows an increase at 
low temperature Jl6[ . The measured ac-susceptibilities 
Xac(T) and Xac{H) exhibit a reentrance both as a func- 
tion of temperature and field The reentrance is ac- 
companied by an out-of-phase response signaling dissi- 
pation and by hysteresis in the field dependence. These 
features are in qualitative agreement with our results for 
the magnetization curve. We find that theory and experi- 
ment agree in order of magnitude for a 2 ~ 0.1, implying a 
transition temperature T C M ~ lOOmK and a spontaneous 
magnetization Mq ~ 1G. A more quantitative compari- 
son with experiment requires a self-consistent treatment 
of the spontaneous currents with the pair potential, ac- 



counting for the nonlocality of the current-field relation 
and its sensitivity to disorder. 

In conclusion, we have demonstrated that the inclusion 
of a finite electron-electron repulsion in a proximity cou- 
pled normal-metal layer naturally produces spontaneous 
interface currents leading to a paramagnetic reentrance in 
the magnetic response: The sign change in the coupling 
across the NS interface leads to the trapping of 7r-states 
at the Fermi energy. The frustrated NS junction relaxes 
through the generation of spontaneous interface currents, 
inducing a paramagnetic moment. A non-trivial issue re- 
mains the requirement that the electron-electron interac- 
tion be repulsive at the low energy scales involved. In- 
teresting consequences of this assumption have been dis- 
cussed in the context of the proximity effect Q and most 
recently in relation to the low temperature transport in 
mesoscopic NS structures 0,0. In fact, the noble metal 
coatings used in the experiments of Mota and co-workers 
appear to be the most plausible candidates for a repul- 
sive electron-electron interaction. Turning the argument 
around, in the light of our findings the experimental ob- 
servation of a paramagnetic reentrance can be taken as 
an indication of the presence of a repulsive interaction. 
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